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It has been shown that the spherically symmetric solutions in a subclass of quasidilaton theory 
are stable against all degrees of freedom and does not even exhibit superluminal propagation. These 
solutions can be found by switching off scalar-tensor interactions, which can not be removed by a 
local transformation. In this paper, we extend the analysis to quasidilaton theory, including non- 
diagonalizable scalar-tensor interactions. We show that all solutions inside the Vainshtein radius 
are problematic : the scalar mode in massive graviton suffers from gradient instabilities, the vector 
mode are infinitely strongly coupled vector perturbations, or the Vainshtein mechanism is absent. 


I. INTRODUCTION 


Modifying gravitational theories typically introduces 
an additional degrees of freedom, which mediates a fifth 
force, and it is tightly constrained by solar system exper¬ 
iments [l[. Therefore, modified gravity requires a mech¬ 
anism to restore general relativity at short distances to 
pass these tests (for a review see [5|). One of the reliable 
mechanism, which reduces to general relativity at short 
distances with an excellent accuracy, could be the Vain¬ 
shtein mechanism [3[, which is originally found in the 
context of the van Dam-Veltman-Zhakarov discontinuity 
3,0 in Fiertz-Pauli theory 0 and can be also found in 
the decoupling limit of Dvali-Gabadadze-Polatti model 

0 and Galileon theory 0. 

Such a mechanism can be seen in ghost-free massive 
gravity called dc Rham-Gabadadze-Tolley (dRGT) mas¬ 
sive gravity 0,13 , whose mass term consists of an in¬ 
finite number of potentials to eliminate the Boulware- 
Deser (BD) ghost 0] (for the proof of the absence of BD 
ghost in full theory, see [0H3)- la the decoupling limit, 
which enable us to capture physics within a certain scale, 
the dynamics of five polarization modes can be described 
within a certain class of scalar(-vector)-tensor theories. 
In the subparameter space such that the scalar-tensor 
interactions can be diagonalized by a local transforma¬ 
tion (this type of theory is called ’’restriced galileon'’ in 
the literature), spherically symmetric assumption allows 
the stable asymptotically cosmological solutions against 
perturbations for all modes although an asymptotically 
Minkowski solution is forbidden due to the appearance of 
hclicity-0 ghost of massive graviton 3- However the in¬ 
stabilities arise when the non-diagonalizable scalar-tensor 
interactions are included in the theory fl8j . 

This fact would motivate us to look for stable solu¬ 
tions in another type of massive gravity. The minimal 
extension of introducing the scalar degree of freedom in 
massive gravity is so-called quasidilaton theory [l9j (for 
a further extension, see 0). This quasidilaton er can 
be introduced by a global symmetry, cf) a —» e a (j) a and 
a —y a — aMpi, where a is a constant and (j> a is the 
Stuckelberg field. We can now ask the same question as 
in dRGT theory : Is there any stable spherically symmet¬ 
ric solutions in quasidilaton theory? Some of the answer 
has been already solved in 3] , which investigated spher¬ 


ically symmetric solutions in the decoupling limit of qua¬ 
sidilaton theory within a certain parameter space where 
non-diagonalizable interactions is absent. The authors 
found that the asymptotically cosmological solutions are 
free of ghost instabilities, gradient instabilities, and su- 
perluminalities although the asymptotically Minkowski 
solution can not be allowed as with the case of restricted 
galileon. Therefore the aim of this paper is to extend the 
analysis in (3 by taking into account the entire param¬ 
eter space of quasidilaton theory. 

This paper is organized as follows. In Sec. II, we de¬ 
scribe the decoupling limit theory of quasidilaton and de¬ 
rive all relevant equations. In Sec III, we find spherically 
symmetric solutions and its consequences to perturba¬ 
tions in the simplest case : no shift-symmetric Horndeski 
terms for a. In Sec IV, we extend Sec III by adding 
shift-symmetric Horndeski terms. Sec V is devoted to 
the summary. 

We adopt the signature (—,+,+,+) for the met¬ 
ric throughout this work, and use the follow¬ 
ing shortcut notation, e ,iapcj E u ^ pa H llv H a B = eHin, 

£i - /OPs j3a UapUfirj = (£?*)£ = 

eeBOA = , and so on. 


II. THE THEORY 


In this paper we consider the decoupling limit of qua¬ 
sidilaton theory. The decoupling limit is defined as 
Mpi —>• 0, m —> 0, A = (Mpim 2 ) 1 / 3 = fixed, and 
T^/Mpi = fixed. All six degrees of freedom (five 
polarization modes in massive graviton and quasidila¬ 
ton) can be decomposed into the massless tensor h^ = 
9nv ~ r W, massless vector A^, and the two scalars 
7 r and a, via the relation <p a = S^x 11 — r] a ^ A^ / (Mp\m) — 
rj^d^-K / (Mp\m 2 ). The scalar-tensor Lagrangian in the 
decoupling limit, up to the relevant energy scale A, can 
be written as 21] 


C = --hreffhap - hr 
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7 1 ? „„„„ ui „ £1 „„ 

+ -f^eennnn - — eeE - -^-eeES 

A 9 12 6 A 3 

+ 2 kr T -' 


(i) 


where is the Einstein operator, II M „ = d^d^n, E pl , = 
dfid^cr, „ is the energy momentum tensor, and a, 05 , 
/3, w, 7 i, and £* are parameters 1 . The Lagrangian C is 
invariant under diffeomorphism transformation h^ —> 
+ 2<9( m £„) and internal Galilean transformations for 


7r and cr, d^i r 


< 9^77 + and d(j7T 


d a n + d,j. The 


terms involving the parameter are obtained from shift- 
symmetric Horndeski Lagrangian [23423 ] for a field 2 , and 
setting all = 0 yields the original quasidilaton theory 
proposed in M ■ If /3 = £4 = 0, which is the case investi¬ 
gated in [U, the action can be recasted into a subclass 
of bi-galileon action [2(il [33 via a local transformation, 
y (cfc/A 3 )7rn^ ly (2£ 2 /A>E^. 

The vector Lagrangian is independent of the tensor 
modes h,, v and quasidilaton < 7 , but couples with the 
scalar 7r [29|, [3(1, 


the equation of motion for yields two independent 
equations for / and a, 


,, 2 M 3 2 

rf = ——-1-A 3 - 2 


M P1 r 


rb! = — 


2 M 

Mpir 


A V 


A - aA 2 - 2/3A 3 
- 2^A 2 - 2£ 4 A 3 
a — (1 + 2aa)X — 6a/3A 2 — 2/3A 3 
- 4&6A. - 66£ 4 A 2 - 2£ 4 A 3 , (5) 


where the prime denotes the derivative with respect to r, 
and we defined the dimensionless variables, 


A = 


A 3 r, 

and the Vainshtein scale, 


A ct = 


\ 47rMpjm 2 / 


A 3 r ’ 


1/3 


( 6 ) 
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where B^ is an auxiliary non-dynamical anti-symmetric 
tensor, and the Lagrangian Ca is invariant under 17(1) 
gauge transformation, A M —> A M + d^x- 

In this paper we consider the general ansatz for spher¬ 
ically symmetric background, 

h 0 o = h(r), hij = f(r)8ij , (3) 

for tensor modes, and 

n(t,x) = ^A 3 f 2 + 7r(r), 

cr(i, as) = ^A 3 t 2 + cr(r), (4) 

for scalar modes. Here we assume that the vector mode 
can be ignored at the background level, A M = 0. Then 


The 7r and a equations of motion can be compactly given 

by 


n-\-m< 5 g 

£ A% m X n X™ = B (i) ( 7 ) , ( 8 ) 

n,m—0 

where * denotes 7r and cr, corresponding to 7r-equation 
and cr-equation respectively, and the coefficients An^ m 
are independent of A and X a , listed in the appendix A. 
On the other hand the coefficients B^ are functions of 
A and A CT , which are given by 

B (7r) = 2(1 + 2aa + 12a/3A + 6/3A 2 ), 

B (CT) = 18£ 4 A CT (26 + A ct ). (9) 

The crucial differences from the case /3 = 0 are A and 
An- dependences on £?b). These terms could potentially 
yield the new type of solutions inside the Vainshtein ra¬ 
dius, and we will see this fact in the next section. On 
the other hand, the asymptotically Minkowski solution, 
which can be obtained by equating the relevant terms 
at the linear regime, is still given by the one in [2l| . 
however 7r is unfortunately the ghost mode around this 
background. Therefore the solutions outside the Vain¬ 
shtein radius should be at least nontrivial cosmological 
solutions, which are A, X a = const although stability con¬ 
ditions need to be investigated. 


1 Note that 7 * are related to the parameters through 70 = (3 — 
4a 5 )/6, 71 = -(2 + a + 2as)/2, 72 = 2(1 + a — as)/3, and 
73 = —(1 + a + a 5 )/6. 

2 The derivation of the decoupling limit of Horndeski theory was 

investigated in )25j . 


III. VAINSHTEIN SOLUTION WITHOUT 
SHIFT-SYMMETRIC HORNDESKI TERMS 


Let us consider the simplest case, £* = 0, correspond¬ 
ing to the absence of shift-symmetric Horndeski terms 
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for a field. In this case, a equation is linear in A CT , thus 
it is analytically solvable for \ a . Then A 5 and A 2 (r*/r ) 3 
in the master equation for 7 r are the dominant compo¬ 
nents well inside the Vainshtein radius, and we have the 
approximate solutions, 

A^zA A^ yi (^y 3 (10) 

r \ r J 

where x\ and y\ are constants, 


X\ 


V i 


± _ fa 

- 2(72 - 4a7 3 ) 2 
2^1(72 ~ 4073) 

LU 



( 11 ) 

( 12 ) 


coefficient A at the leading order inside the Vainshtein 
radius, 


3 c 2 n 

1 / 2(72 - 4 a 7 3 ) 2 + P 2 uj 
32 u V Pl3 



A (17) 


Since u > 0, one cannot eliminate the leading order of c 2 
by any combination of the parameters. Thus we arrive at 
gradient instabilities in both radial and angular direction 
for 7r perturbations. 

Another type of solution can be obtained by setting 
the right-hand side in Eq. © to be zero, = 0, or 
explicitly, 


1 + 2aa + 12a/?A + 6 /?A 2 = 0. (18) 


Then the equations of motion for / and h can be rewrited 
as 


(1 + 

A" + °( 

})• 

(13) 

■(1 + 

0 . 2 M , 

n\ 


Xlfj) M P ,r 2 ' 

a) 

(14) 


One can clearly see that the contribution from helicity-0 
mode cannot be screened in this solution, i.e., there is 
no Vainshtein mechanism. However we can redefine the 
Plank mass (or Newton’s constant) as the one that we 
observe at short distances such as M P \ = Mpi/(1 + x\(3) 
so that the leading order in the gravitational potential 
agrees with Newtonian one. However, the parametrized 
post-Newtonian expansion gives the different result as 
one in general relativity, which could be tightly con¬ 
strained by solar system experiments. Therefore, we 
require x\ft <C 1. Now let’s take a look at the scalar 
perturbations around this background. The quadratic 
Lagrangian for the scalar perturbations are given by 

£ (2) D A(d t $) 2 + B(dt^) 2 + C(c^)(dA0, 

= A^d t ij) + + (b-^ j(d t ip) 2 (15) 

where <p(t,x) and ip(t,x) are the perturbations of 7 r and 
a, respectively, and the coefficients A and B are given by 

A ~ -127 3 x 2 yi (16) 

and C oc (r*/r) , making C 2 /AA the small correction in 
Eq. (fl5l) . Thus we require —7 3 ^f J/i > 0 to avoid ghost 
for (diagonalized) 7 r field and w > 0 for a field. This 
feature is completely contrast to the case /3 = 0, which 
is C 2 /4A B inside the Vainshtein radius, leading to 
the ghost mode for the scalar perturbations. Next we 
would like to evaluate the sound speed of 7 r perturbations. 
Following to the method in [2l|, one can find that the 
radial and angular sound speeds are related to the kinetic 


In this case, A is constant everywhere, and the 7 r force 
is successfully screened since the additional contribution 
from the scalar is S(f'), S(h') ~ const . Then the leading 
contribution to the gradient energy in 7r’s perturbation 
is given by 


£(2) 


—3 (3(a + A) 



3 


2(d r (j)) 2 


(dot) 2 


(19) 


If A ^ —a, one of the squared sound speeds of radial 
and angular part is always negative, leading to the gra¬ 
dient instability in the 7 r-sector. This feature due to non- 
diagonalizable interactions is already addressed in the 
context of the decoupling limit of dRGT massive grav¬ 
ity 0 and Horndeski theory [25j. Therefore we require 
A = —a at short distances. However the kinetic term of 
the vector perturbations are given by 


p(2) 

l 'A 


1 - 2aA - 6/3A 2 
4 -(- 2 a — 2A 


(9* A) 2 + • • •, 


( 20 ) 


where we set the gauge condition V • A = 0 and A M = 
(0, A) 3 . One can read off the numerator of the coefficient 
is the same combination that we imposed in (1181) if A = 
—a. Thus the vector perturbations is infinitely strongly 
coupled in this solution. 

It should be noted that the scalar field 7 r is the Lorentz 
invariant form 7 r = (a/ 2)x IJ, x fl if A = —a, and in this 
case the stable self-accelerating solution can be found in 
a broad parameter range [3l[. However, the condition 
m is not the case of the stable one because the vec¬ 
tor perturbations are infinitely strongly coupled. If the 
solution asymptotically approaches this stable de Sitter 
space-time, the solution inside the Vainshtein radius is 
described by (flOl) and the 7 r perturbation suffers from the 
gradient instabilities discussed in the above. Therefore 
the self-accelerating solution found in [UJ is problematic 
at short distances. 


3 The auxiliary tensor B /n y is already integrated out in this ex¬ 
pression. 
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IV. VAINSHTEIN SOLUTION WITH 
SHIFT-SYMMETRIC HORNDESKI TERMS 

In the case /3 = 0, the shift-symmetric Horndeski terms 
for a field are crucial for the stable perturbations and 
subluminal propagations [21]. In this section we include 
the shift-symmetric Horndeski terms in the case f3 ^ 0. 
In the presence of these terms, i.e., ^ 0, the equations 

for 7r and cr are coupled quintic equations, thus we assume 
the following ansatz inside the Vainshtein radius, 

A~xi —, ~yi — . (21) 

r r 

Then the equations of motion for 7r and a gives the same 
equations for x\ and y\, 

1 + /?£?+£ 4 ^ = 0 . ( 22 ) 

By using these solutions, the metric fluctuations can be 
written as 

/' = -(1 + f)xl + Uv\)-^s + o(i), 
h’ = + M + + (23) 


V. SUMMARY 

In this paper we investigated the possibility of stable 
spherically symmetric solutions in the decoupling limit 
of quasidilaton theory in the whole parameter space. We 
showed that the presence of non-diagonalizable scalar- 
tensor interactions contains the following serious prob¬ 
lems. One of the solutions inside the Vainshtein radius 
can not be allowed due to the appearance of the gradi¬ 
ent instabilities for 7r perturbations. For the solution, 
which does not depend on the source term in equations 
of motions, the extra degrees of freedom 7r and cr can be 
successfully screened, but the vector perturbations are 
infinitely strongly coupled for any parameter space. We 
confirmed that these conclusions can not be evaded in 
the case of inclusion of shift-symmetric Horndeski inter¬ 
actions for <7 field. One of the solutions does not even 
have the Newtonian gravitational potential at the leading 
order due to the cancellation with the contributions from 
the scalar modes. The other solution, whose equations of 
motion is independent of the source term, encounters the 
same problem as in the case of the absence of Horndeski 
terms. Therefore, the case f3 = 0 found in [21] is the only 
consistent quasidilaton theory, which are free of ghosts, 
tachyons, gradient instability, and superluminality and is 
not ruled out by solar system experiments. 


One can clearly see that the leading term are canceled 
due to Eq. C^l) , and both scalar modes completely screen 
the term from massless graviton. Then the leading terms 
in these metric fluctuations are f ~ h! ~ 1/r; therefore, 
this solution can not even reproduce the Newtonian pro¬ 
file. We disregard this solution for this obvious reason. 

As in the previous section, we have the other type of 
solution, A, A ct = constant everywhere in space. For 
(i / 0 case, we further impose the condition B^A = 0 
in addition to = 0, which completely eliminates the 
source terms in both ir and cr equations. Then these 
conditions translate into Eq. m and 

A CT (26 + A ct ) = 0. (24) 

As one can see, introducing the shift-symmetric Horn¬ 
deski does not change anything about the condition (fTSl) 
and the gradient energy of 7r’s perturbations CHI as well 
as the vector perturbations (because the vector mode 
only couples with 7r). Therefore we conclude that the 
vector perturbations are infinitely strongly coupled even 
in the presence of Horndeski terms. Furthermore, the 
leading contribution to the gradient energy for cr’s per¬ 
turbation is given by 


4 2) = -3^4(6+ 


3 r 


2 (drip ) 2 - (< 9 fiV 0 2 


(25) 


We have X a = 0, —26 from Eq. (liMl) . which means that the 
<7 perturbations always suffer from gradient instabilities. 
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Appendix A: Coefficients in equations of motion for 
7r and a 

In this appendix we summarize the coefficients of equa¬ 
tions of motion defined in ([ 5 ]). A^m are given by 

4*o = -a+ 467 o, 

= 3 + 6 acx + 8671, 

A^l = — 6 (a + aa 2 — 4 a/3 — 2672 ), 

A^g = 2(a 2 — 4/3 — 20aa/3 + 8673), 

4g = —60a/3 2 , 

4To = - 12 / 3 2 , 

Agj = - 4 ( 37 0 - 2071 + 6£ 2 ), 

A$ = 2(26+200^2-3^4), 

4 $ = 2(1 + 2 aa)U, 

Aqi = -8(271 - 3072 - 6a£ 2 )> 
a [*2 = - 4(«6 — 6 a / 3£ 2 — 36 a ^ 4 ), 
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41 = -24a/3£ 4 , 

41 = -12(72 - 4a73 - 26/36), 

Ag = — 36&/3^4, 

41 = —12/36, (Al) 

and A^i are given by 

A <1 = 4a 5 + 6070 - buj, 

A$ = -6(370 - 2 a 7 i - 66), 

A 2 4] = -3(471-6072 + 6^2), 

41 = -6(72-4073 + 26/3^2), 

A 07 = 3(w-2&6+2a6), 

41 = 3(2^-126^-66^ + 30^4), 


41 = 6(2^+6- 2066^4-26^), 

41 = 906^, 

4:1 = ~ i8 & 

Ail = ~ 6 ( 2 & + aa ^2 - 366), 

A^l = —3(3 + 2aa)6, 

4:1 = 3(a6 - 12a/36 - 26a6), 

41 = -54a/36, 

41 = -366/36, 

41 = -18)86, (A2) 

and any other coefficients are zero. 
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